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Abstract 

We analyze a simultaneous continuous measurement of photon counting and homodyne detection. 
The stochastic master equation or stochastic Schrodinger equation describing the measurement 
process includes both jump-type and diffusive- type stochastic increments. Analytic expressions 
of the wave function conditioned on homodyne and photocount records are obtained, yielding 
the probability density distributions and generating functions of the measurement records. The 
obtained results are applied to typical initial conditions — coherent, number, thermal, and squeezed 
states. Monte Carlo simulations of the measurement processes are also presented. 
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I. INTRODUCTION 



Since the birth of quantum theory, the problem of quantum measurement has been among 
the most thought-provoking problems in physics. In the early stage of quantum theory, the 
postulate of projective measurements, also known as von Neumann measurements, was put 
forth by von Neumann [1]. The von Neumann's postulate, however, is not applicable to most 
experimental setups in a straightforward manner. A more general theory of quantum mea- 
surement which can describe non-pro jective measurements such as photon-counting measure- 
ment [2] was formulated much later [3], and it was clarified that the theory of measurement 
is closely related to the ideas of open quantum systems and quantum operations [4, 5]. 

Many open quantum systems are known to obey a special form of master equation, called 
the Lindblad master equation In a seminal paper [6], Lindblad showed that this master 
equation can be derived from the mathematical requirement that the quantum operational 
map take the semigroup form exp(t£). From the view point of continuous measurement, 
this time evlolution for the density operator corresponds to a situation in which we discard 
the measurement outcomes, (see Ref. [7] for a review). 

The theory of the continuous quantum measurement has been studied from various stand- 
points. One is a mathemaical approch [8-10]. There are also experimentally oriented pieces 
of work especially in quantum-optical systems [2, 14-17] and, more recently, in mesoscopic 
systems [18]. 

There are two distinct types of measurement processes: a jump-type and diffusive-type 
processes. In the jump-type measurement, discontinuous state changes occur at discrete 
times. A typical example is the photon-counting measurement, where the coupling between 
the photon field and the detector is adjusted so that the probability of more than one 
photon being detected during any infinitesimal time interval (i.e., the resolution time) is 
negligible. Therefore, the photon-counting measurement consists of two fundamental pro- 
cesses: no-count and one-count processes. The state change in the no-count process is, 
however, different from that of measurement-free evolution due to the back action of the 
measurement. On the other hand, in the diffusive- type measurement, the state change is 
continuous. Examples include balanced homodyne measurement [15] and continuous obser- 
vations of the position of a particle [8]. In classical stochastic processes, the jump-type and 
diffusive-type continuous measurement correspond to the Poissonian and Wiener processes, 
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respectively [19, 20] 

Mathematically, the hybrid type of continuous measurement with diffusive and jump out- 
comes is also possible. A general equation of the continuously observed system was derived 
under general semigroup assumptions by Barchielli et al. [10]. The simultaneous measure- 
ment of photon counting and homodyne detection analyzed in this paper is an example of the 
hybrid type continuous measurement. We note that a simultaneous measurement of photon 
counting and homodyne detection was discussed in Ref. [11, 12] in a different context. We 
also note that continuous measurements driven by Levy processes are discussed in Ref. [13] 

This paper is organized as follows. In Sec. II, we apply the general theory tothe simul- 
taneous measurement of photon counting and homodyne detection. We derive an analytic 
expression of the conditional wave function. In Sec. Ill, we examine probability laws of 
measurement records by deriving the probability distributions and the generating functional 
for measurement outcomes. In Sec. IV, we apply obtained general expressions to typical 
initial quantum states, namely, the coherent, number, thermal, and squeezed states. We 
also present the results of the Monte Carlo simulations for each of these initial states to 
illustrate how the hybrid-type measurement changes the average photon number. In Sec. 
V, we summarize the main results of this paper. In the Appendix, we show derivations of 
some formulae used in the main text. 

II. HOMODYNE MEASUREMENT ACCOMPANIED BY PHOTON COUNTING 

In this section, we consider a simultaneous measurement of photon counting and homo- 
dyne detection. 

A. Setup of the system 

The measurement scheme discussed in this section consists of photon counting and bal- 
anced homodyne detection. A single-mode photon field, described by the annihilation oper- 
ator a, is divided by a beam splitter into two, one of which is detected by a photodetector 
and the other is superimposed by a local oscillator with amplitude and then measured by 
a balanced homodyne detector, as schematically illustrated in Fig. 1. 

The measurement process during an infinitesimal time interval dt can be mathematically 
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vacuum field b 




FIG. 1: Schematic illustration of the system. An input photon field is simultaneously measured by 
a photodetector and a balanced homodyne detector. See the text for details. 

expressed by following measurement operators: 



where 71 > and 72 > denote the coupling strength of the photodetector and that of 
the homodyne detector, respectively. The unitary part of the time evolution is given by 
H = wn, where u is the detuning of the photon field with respect to the local oscillator /3 
and n := at a is the number operator of the photon field. The Planck constant h is set to 
be 1 throughout this paper. dW is the stochastic variable corresponding to the homodyne 
record with the expression [7] 



where (■) := tr[p-] and dW is the Wiener increment which obeys the Ito rule (dW) 2 = 
dt [13, 21, 22]. Equation (1) describes homodyne detection, while the measurement opera- 
tor in Eq. (2) corresponds to the photodetection event. Actually, the vacuum field enters 
the mirror before the photodetector, but this effect can be neglected because it does not 
contribute to photodetector. 




(1) 



Mi = a/ jidta, 



(2) 



dW = (a + d ] )dt + dW, 



(3) 
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B. Wave function for the no-count process 



Let the initial state vector of the system at t — to be l^o)- We consider the time 
evolution of the wave function under the condition that the homodyne records are W(-) and 
that there is no photocount. It is described by the measurement operator in Eq. (1). Thus, 
the unnormalized wave function | during the no-count process obeys 



m+dt)) = M dW m)) 



I iu) + — ) hdt + ^J^MW 



(4) 



where T := 71 + 72 and the tilde over ip indicates that the state vector is unnormalized. To 
solve this equation, let us introduce an interaction-picture wave function 



\ipi(t)) := exp 



(t-t )[iu + -]n 



\m). 



(5) 



From Eq. (4), the time evolution of \ifji(t)) is given by 

d\Mt)) = ^e-i^y^ddWtlMt)). 



(6) 



Therefore, we obtain 



\ipi(t + dt)) = exp 



72 e 



d^&dWt - ^ e -(«"+r)(*-*o) a 2 dt 



!&(*)>■ (7) 



Note that the second term in the exponent arises from the Ito rule dW 2 = dt. Since the two 
operators in the exponent of Eq. (7) commute with each other at any differnt time, Eq. (7) 
can be solved iteratively, giving 



\i>i(t)) = exp 



exp 



72 / e 

t 



K|)(t'-to)^ a _ 72 I* e _ (2iw+ r)(f-to) df / a 2 

2 J t 



1 v\, , s ~ 1 - P -(2^+r)(t-t ) 



to 



2iu + T 



-a 



l^o) 
l^o). 



Going back to the original \ip(t)), we obtain the unnormalized wave function for the no-count 
process: 

-t . j 2 1 _ e -(2tw+r)(t-t ) 



\$(t)) = e -H^) ( *-* o) "exp 
=: N(t,t ;W(-))\i> ). 



72 / e 

to 



-(*-+|)(f-to) tfl y a _ 



2iw + T 



-a 



l^o) 



(9) 
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C. Wave function for the m-count process 



We generalize the result obtained in the previous subsection to the conditional wave 
function under the m-count process. Let us assume that the initial condition of the wave 
function is \ip(t = 0)) = |^o)- The conditional wave function under the condition that 
photocounts occur at times t±, t 2l ■ ■ ■ ,t m during the time interval (0, t) (0 < t± < t 2 < ■ ■ ■ < 
t m < t) and that the homodyne records are W(-) is given by 

m)} = (dtr^m-, t u t 2l ... t m] ww), m 

where 

t 2 ,...t m ;W(-))) 

= ^ n/2 N(t, t m , W(-))aN(t m , t m _ x - W(-))a ■ ■ ■ aN(t 2 , h; W(-))aN(t u 0; W(-))\tf> ). (11) 
The product of the operators on the rhs of Eq. (11) can be simplified as 

N(t, t m , W(-))aN(t m , t m _ i; W(-))a ■ ■ ■ aN(t 2 , h; W(-))aN(t u 0; W(-)) 

= e -(^+|)(ti+t 2 +-+t m ) e -(ia,+E)tft a m exp [ At& + 5 ( t ) fi 2j ? ( 12 ) 

where A t and B(t) are given by 

A t :=V72 fe-^WdW*, (13) 
J o 

Bit) := -Jl— . (14) 

Thus, the conditional m-count wave function is given by 
\4>(t;h, t 2 ,...t m ;W(-))) 

= 7l m /2 e -(- + £)( tl+t2+ ... +tm ) e -(^+E)tn amexp ^ + B{t y a 2^ ^ 

Note that the time dependence on photocounts of this wave function arises, aside from the 
c-number factor e ~( itJ+ 2)( <1+ * 2+ '+' m ) ; only through the number of photocounts m. 

D. Time Development of Expectation Values of Observables 

For later use, we derive the equation for conditional time development of expectation 
value (A) of an arbitrary operator A. 
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For a one-count process, the density matrix p(t + ) immediately after the state p(t) at 
time t is given by p(t + ) = apa)/(n). Thus, the expectation value (A) + immediately after 
the photocount event is given by 



(A)- 



(h) 



(16) 



For a no-count process, the conditional evolution of the density matrix with homodyne 
record dW is given by 



t 



p + dp 



M awPM dW 



(17) 



Substituting the expression of the measurement operator in Eq. (1), we obtain the equation 
for the differential of the expectation value (A) in our system as follows: 

tr[AM dW pMl w ] 



d(A) 



tT[M dW pM\ 



{A) 



dW 1 



MM) - y({Ai, Ah}) - f({A,h} - 2aUa) 
+ ^(AiAa + AtfAA)dW, 



72 



dt 



(18) 



where A A := A - (A). 

Important examples of the expectation value (A) are average photon number (n) and 
quadrature amplitude (a). By substituting h and a into A, we obtain for the one-count 
process 



(n>. 



(n 2 ) — (h) 

W) 



(n> 



([An] 2 ) 
(n) ' 



. . (a) dd) (Ah Ad) 

' a )+ = ^TZV 1 = (<>>) + 



(h) 



(h) 



and for the no-count process 



d{h) 
d(d) 



{j2{n) + 7i([Ari] 2 )) dt + ^{AhAd + Ad ] Ah) dW, 



dt + ^({Ad} 2 + Ad ] Ad)dW. 



%uj + - ) (a) + (AnAa) 



(19) 
(20) 

(21) 
(22) 



III. PROBABILITY LAWS OF MEASUREMENT RECORDS 



A. Probability density functions 



In this section, we will derive general results on the probability distributions of homodyne 
and photocount records. From the general considerations on the measurement operators, the 
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joint probability density of homodyne records W(-) and photodetection times ti,t 2 , ■ ■ ■ ,t m 
is given by 

dhdh-'-dtm x (io(W(-)) x $(t;ti, h,...t m -W{-))\^{t-t 1 , t 2 ,...t m ;W (■))), (23) 

where /j,q is the Wiener measure. The square of the norm of the wave function is evaluated 
to be 

t 2 ,...t m ;W(-))\iP(t;t 1: t 2 , . . . t m , W(-))) 
= 7l m e- r(tl+t2+ - + * m) (^ | exp [Ald ] + B*(t)(a t ) 2 ] (a t ) m e- rt "a m exp [A t a + B{t)a 2 ] |Vo) 

_ 7 m e -r(ti+t 2 +-+t m )^ o j . e ylia+A*at+B(i)a 2 +B*(i)(at)2_(l_ e -ri)ata^ t ^ m . ^ 

where the symbol : • • • : in Eq. (24) indicates normal ordering which places annihilation 
operators to the right of creation operators. In deriving the last equality in Eq. (24), we 
used the formula 

e xh = . e ata( e --i) . ^ 25 ^ 

which is valid for an arbitrary c-number x. In the limit of t — > oo, Eq. (24) reduces to 

T m e -r(t 1 +t a +...+t m )^| . e A 0o a+Ai,at + B( 0O )a3 + B.( 0O )(at)3_ata^t^ m . |^ ^ 

where 

poo 

A 00 = ^T 2 e-V"+ T *YdW v , (27) 



o 



B(oo) = (28) 

v ; 2zw + r v ; 

The joint probability p m (t; W(-)) of m-photocounts being recorded during time interval 
(0,t) and the homodyne records W(-) is obtained by integrating Eq. (24) with respect to 
h, t 2 , ■ ■ ■ ,t m in the integration range < t\ < t 2 < • ■ ■ < t m < t. The relevant part of the 
integral is the exponential e~ r (* 1+ * 2+ '"' m ) and evaluated to give 

ft ft m ft2 

/ dt m / dt m ^ ■ ■ dhe-^^-^ 

Jo Jo Jo 
1 /i_ e -rt\™ 

= ;sHH ■ (29) 
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Thus, p m (t; W(-)) is given by 
Pm(t;W(-)) 

= U Q \ : J_ f2l(i _ e -r i)a t a y m e A t a + A*at +B(i )a 2+ B. W (at)^(i- e --)ata . ^ (3Q) 
m! VT / 

3 A t *at+B*(i)(at)2 [ n \ p -ri(n-m) fTl/-, _ -Ft s\ m A t a+B{t)a* ■ 



(V>o|e ' ( n } I ^-^"-^(^(l-e- 1 *)) e A ' a+ ^ a >o), (31) 
mj 

where in deriving the last equality we have used the formula 

(oJ) m e xfl a m = h(h - 1) • • • (n - m + l)e x{h - m \ (32) 
and defined the operator binomial coefficient by 

n | h(n — 1) • • ■ (n — m + 1) 



mi 



(33) 



m 

In the limit t —> oo, the joint probability in Eq. (30) reduces to 

<Vo] : A fS^aTexp |X,a + A*^ + 5(oo)a 2 + fi*(oo)(a t ) 2 - a)a] : |^>- (34) 
ml V 1 / L J 

Note that Eq. (26) gives the total probability functional for the measurement outcomes 

and that its dependence on the homodyne records enters this formula only through A^. 



1. Marginal distribution of homodyne records 

The marginal (unconditional) distribution of homodyne records p(t; W(-)) can be ob- 
tained by taking the sum of p m {t] W(-)) over the number of photocounts m. This calculation 
precedes as follows: 

P{t-W{-)) 

oo 
m=0 

= U Q \ : J_ (llh - e -r i)a t a r e A t a + A*at +B(i )a 2+ ^ W (at)^(i- e -")ata . , - 

ml \ r / 
= (^1 : e Ata+^at + B(t)a=»+B-(t)(at) a _^ ( i_ e -rt )a t a . ^ (35) 

which, in the limit of t — > oo, reduces to 

p(oo; W(-)) = : e ^a+^at+ij(oo)a2 + B*(oo)(at)^^ata . ^ 

This functional of W(-) gives the marginal distribution of the homodyne records. 
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2. Marginal distribution of photon- counting records 



We can also derive the marginal photocount distribution of p m (t), which is the probability 
of m photocounts being registered in the time interval (0, t) and can be derived by averaging 
p m (t; W(-)) over Wiener measure fj, (W(-)). To evaluate the integral, we use the following 
general formula on the integration of a normally-ordered operator-valued function: 

f(a,a j ; u) : fi(du) =: J f (a, a j ; u) fx(du) :, (37) 

where /(a, a*]u) is a function of operators, a and a', and of integration variable u. The 
integration on the rhs of Eq. (37) can be carried out as if a and a) were c-numbers. By 
using the normally ordered expression of p m (t; W(-)) in Eq. (30) and the formula (37), the 
relevant part is the linear term on A t : 



Ata+Afai 



e - ?*' { e - iut ' &+e iut> tf)dW t t 



By using the formula 



the integral on the rhs in Eq. (38) becomes 



exp 



72 ( 1 - e-^+n* ^ 2 1 - e -(-M"+r)* ^ | 72 



2 V 2iu + T 
Thus, the photocount probability distribution is given by 

1 /7i 



A / 7i \ 
Pm(t) = (V'ol : — r -e" r *)a f a exp 

m! V 1 / 

One may use Eq. (25) to rewrite Eq. (41) as follows: 



r v ' 



■■ IV'o). 




Pm{t) = <*0o| 

In the limit of t — > oo, Eq. (42) reduces to 

p m (oo) = (^0 

= (^0 



-rtN 



7iC 



-rt 



72 



l^o). 



^(F ataJ 



7i. + 
— a T a 



l^o) 



(r) l^>- 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 
(44) 
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B. Generating functional 



In this subsection, we derive a general formula for the generating functional of measure 
ment records dW t and dN t , where dN t is defined by 

otherwise. 



dN f 



(45) 



Instead of deriving the generating functional of dW t and dNt, we discuss that of dA t 



7 2 e y 2+lu) rdW t and dN t , which is defined as 



e f °° t{t')dA t ,+!™ t*{t')dA* t , +/ °° v (t')dN t , 



(46) 



where £ and t] are arbitrary functions of t. This functional contains all the information 
about the probability distribution of measurement records. 

To calculate this generating functional, we note that the stochastic integral J °° rj(t')dN t i 
becomes ^2^=1 vi^k), if the photocounts occur at times t x < t 2 < ■ • • < t m . Thus, the 
generating functional can be evaluated as 

m=0 Jo Jo Jo J 

x (Jaofah,--- ,t m ;W(-))|^oo(ti,t 2 ,--- ,t m ;WO)}e^W^U^^(0^ + Er=o^). (47) 
From the square norm of the wave function in Eq. (26), we obtain 

Af[«o,r(o,i(-)] 



£ 

m=0 



(II 

JO 



dt 



m—l 



*2 



dh / //o(W(-))7i m e 



"i„-r(ti+i 2 +-+t m )+Er=o T ?(*fe) 



771=0 ' ^ 



fl tfa J e- rt ' + ^dt' 



x e V^ I °° dW t ,[e K^+V* «(f)+a)+e I (r(t')+a t )]+-B(oo)a2+B*(oo)(at)2-ata . 



exp 



- / ie(f)e~ (iw+ ^ t '+r(Oe~(~ iw+ ^ l/ rrft' 

2 Jo 



(V> | : e 



where 



« = 72 / le 

'0 

r>00 



z/ = 71 



(2-+r)^ (t , ) + e -rt' r(f))( ^ 
l)<ft'. 



e -r*' (e r,(*') 



(49) 
(50) 
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Equation (48) gives the general formula for the generating functional of measurement 
records. 

The generating function E[e^ At+ ^ A ' +r,Nt ] with respect to output variables A t and N t can 
be derived by substituting 

m=\ (51) 
I (otherwise), 

V (t>) = I' (52) 
I (otherwise) 

into Eq. (48). The result takes a form similar to Eq. (48): 

fil^At+Z* A* t +r,N t ^ 
72 fl ~ 



exp 



1 _ p -{2iuj+T)t' i _ -(-2iw+T)t i _ -rt 

e + . , . (D 2 + 2^ — — 



2 V 2iuj + F -2icu + T 

x (^ol : exp[n t a + n* t a) + z/ 4 a^a] : |^ ); (53) 



where 



e -(2iuj+T)t' l - e 



rt 



"■ = n 2,wr ^ + (54) 

„, = 2l(l- e - r ')( e '>-l). (55) 

To gain the physical insights into the formulas (48) and (53), let us assume that a and a) 
were c-numbers. Then, the generating functional (48) would be that of independent stochas- 
tic processes, dW t and dN t . Such a description is justified only when the initial state is a 
coherent state (see Sec. IV); otherwise, there will, in general, be correlations between these 
output records because of the noncommutativity of a and a) . In this sense, the measure- 
ment records give us information about the system's deviation from the coherent state as 
exemplified in the next section. From a stand point of measurement theory, the correlation 
between the two output records, which can be seen in the generating functional (48), reflect 
the backaction of one measurement channel on the other. 



IV. APPLICATION TO TYPICAL INITIAL CONDITIONS 



In this section we apply the general formulae obtained in the previous sections to typical 
quantum states: coherent, number, thermal, and squeezed states. 
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A. Coherent state 



The coherent state \a) is represented in the number-state basis as 

, , _mL or . . . . 

\ a ) = e 2 Z^^tI 71 )' ( 56 ) 



where a is an arbitrary complex number. The coherent state is an eigenstate of the boson 
anihilation operator: 

a\a) = a\a). (57) 

The conditional wave function of the m-count process in Eq. (15) for an initial state 
l^o) = \ a ) is given by 

h, t 2 , . . . t m , W(-))) =7 -/2 e -(-+i)( il+i2 +---+^) a m e A t a + B(^-^(l- e ^ t) 

x {e-^+^a). (58) 
In evaluating the m-count wave function, we have used Eq. (57) and the formula 

e A >> = e-^ 1 "^ V«>- (59) 



It follows from Eq. (58) that the normalized state vector does not depend on the measurement 
outcomes and is given by \e~^ + ^ tn a); the system develops deterministically. 
The generating functional in Eq. (48) is evaluated as 

M[£(-U*(-),^(-)] 



poo 

— / uty-^y + cit')e<- iw+ ^) t, \ 2 dt' 

2 Jo 



exp[na + K*a* + u\a\ ], (60) 



= exp 

where we used the formula 

(a I : /(a,a f ) : |a) = /(a, a*). (61) 

Note that the generating functional in Eq. (60) can be obtained by replacing a and a) by 
the corresponding c-numbers, a and a*. From the definitions of k and v in Eqs. (49) and 
(50), we find from Eq. (60) that 

MU-),C(-)M-)\ =M[£(-u*(-),^(-) = o] 

xM[£(-)= 0,770], (62) 
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which implies that homodyne records W(-) and photocount records are statistically inde- 
pendent. 

From Eq. (53), the generating function of A t is 



"72 (I ~ e (2lbJ+T V 



exp 



2iu + T 



1 _ p -(-2iuj+T)t 1 _ -rt 

e 2 + — ~ — —(cf + 2^- e 



-2iu + r 



+£72 a 



I _ e (2iw+r)t' 



1 * 
+ a 



1-e 



-rt 



T72 a 



1-e 



■lei- 



(-2iaj+r)t' 



+ a- 



1-e 



-rt 



(63) 



2iw + T T / V -2iu + T 

which implies that A t is a complex Gaussian variable with its first and second moments 
given as follows: 



E[At] 


( 1 

= 72 ( OL 


E[{A t 


-E[A t ]f] 


E[\A t 


-E[A t f] 



f a* 



rt 



2iu + r r 

72 1 _ e (2^+r)t' 

~ ~2 2iu + T ' 
1 - e~ rt 



(64) 

(65) 
(66) 



From the normally-ordered representation (41), the marginal photocount distribution is 
evaluated as 



ml \ 1 



m 

Q!|" I exp 



^(l-e- r ')l«| 2 



which is a Poisson distribution with its mean ^-(1 — e r *)|a| 2 . 



B. Number state 



(67) 



The m-count wave function for an initial number state |^o) = \n) is evaluated as 

t 2 ,...t m ;W(-))) 

= T -/ 2 e -(-+0(*i+ta+-+^) e -(-+l)"*exp [A t a + 5(t)a 2 ] a m \n) 

% Lfc/2J ^- 2 '5(t)' 



^W2 e -(^+E)(ti+t 2 +-+t m ) e -(iaj+|)nt y2 



k=0 1=0 



(k-2l)W. V (n-m-k)\ 



\n — m — k) 



^W2 e -(ta;+5)(ti+t2+-+tm) 



-m L fc / 2 J Ak-2l-n/.\l 



X 



EE 

fc=0 z=o 



(/c - 2/)!/! V {n-m-k)\' 



-(iw+E)(n-m-fc)t| n _ m _fc^ 



(68) 



14 



where [^J is the largest integer that does not exceed x. Note that m < n since 
\ip(t; t\, t2, ■ ■ ■ t m ; W(-))) vanishes for m > n. 

The joint distribution function p m (t; W(-)) in Eq. (24) becomes 



Pm(t;W(-)) 



(1-e 



(1-e 



-Ft- 



-Ft- 



I — {iui+-^)tfi 



E 

k=0 



exp[A t a + B(t)a 2 ]a m \n) 

2 



Lfc/2J At 2l B( t y 



5 ( fc - 2/ )!^ 



(n — m — k)\ 



-Ft(n—m~k) 



(69) 



In the limit of t — > oo, Eq. (69) reduces to 
p m (oo; #"(•)) = n! 



^ (ra-m-2Z)M 



1=0 



(70) 



n! 



The generating functional takes a simpler form. The calculation preceeds as follows: from 

n 

e ™\n) = 

we have 



m=0 



m\ y (n — m!) 



|n — m), 



(n| : exp[fi;d + K*d^ + z/d^d] : |n) 
= {e K& \n))\\ + i/)V a |n) 



E 

m=0 



'm!) 2 (n — m)! 



< m (l + z/) 



(1 + v) n L n 



AC 



1 + V 

where L n (x) is the Laguerre polynomial defined by 



") = E 



nl 



n\dx ny ^-^ (ml) 2 (n — m)\ 

m=0 K ' v ' 



Thus, the generating functional and the generating function are given by 
M[£(-),f(-)^(-)] 



= exp 
= exp 



/•LXJ 

- / l^'K^^' + e(t')e<-^y\ 2 dt' 
2 Jo 



(l + v) n L n 



1 + v 



ry n f\ _ p -(2iw+r)t' i _ -(-2iw+r)t i _ -rt 



x (1 + v t ) n L r 



(71) 



(72) 



(73) 



(74) 



(75) 
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Note that for the initial number state the generating function(al) cannot be factorized into 
a function of £ and that of 77, implying that there are correlations between these two mea- 
surement outputs. 

From Eq. (42), the marginal counting distribution is 



which is a binomial distribution. 

The Monte Carlo simulation for the measurement process is done for the number-state 
initial condition. The simulation method is as follows: For each time step At, we first check 
if the photodetection occurs or not. If it does, the state vector or density operator evolves 
according to the jump operator in Eq. (2). If not, the state evolves according to the diffusive 
measurement operator in Eq. (1) corresponding to homodyne detection. 

The results for an initial number state \n) with n = 3 are shown in Fig. 2 (a). In the 
no-count event, the expectation value of the photon number decreases in average according 
to Eq. (21), while there are local stochastic deviations which arise from the diffusive term 
^(AnAa + Aa) Ah)dW in Eq. (21). 

C. Thermal state 

Let us assume now that the initial state is a thermal state 




(76) 



—/3ujh 



00 
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(l-e^)^e 



—fium 



n)(n 



n=0 



where 



Z = tr[e-^™] 



1 



1 - e-P" ' 
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and (3 is the inverse temperature. Then, p m (t; W(-)) can be calculated by taking the ensemble 
average of the corresponding quantity for the initially number state over n: 



oo 

45> 



(— (1 - e" rt ' 



E 

fc=0 



^ (*-2Q!l! 



-Tt(n— m— k) 



E 

fc=0 



7^ ( k ~ 2l ) m 



(n — m — fc)! 



(77) 



where in the last equality the following formula was used: 

E(n + m)\ n x k (k + m)\ 



n=0 



(n-k)\ (l-x) m+k+1 ' 



(78) 



The generating functional and the generating function are evaluated as follows: 



M[£(-),f(-)^(0] 



cxp 



72 |^_ t/ ^_('i w+ ry ( t *^„-f-iw+^t'|2^ 
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1 e -/3ojn 



1 + ^) n L r 



n=0 



1 + 1/ 
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- 1 - v 
- 1 

" epu - i _ Vt 



exp 



72 / " i^-L+rii- 
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+ ■£ ie(Oe" 



, (79) 



+ 



72 



1 _ e -(2iw+r)f ]_ _ e -(-2iw+r)< 
< 2 + - 



*\2 



1 - e 



-Ff 



■ler 



2 v 2iw + r -2iu + r y ' r 

In deriving the last equality in Eq. (79) we used the relation 

°^ i \ xt ' 

2_^t n L n (x) = j—^exp ~Y~t 

n=0 L 

Again, these characteristic functions are not separable with respect to £ and r\. 
The marginal photocount distribution is 



(80) 



(81) 



Pm(t) = 



-fiuin 



n 



p(t) m (l-p(t)) n 
Pit) 

- 1 + p(t) -1+ p(t) 



n=0 



(82) 
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(a) (b) (c) 




FIG. 2: (Color) Monte Carlo paths of (n) starting from (a) number, (b) thermal and (c) squeezed 
states. In each figure, the dashed curve shows the ensemble average of the photon number which is 
given by (h)oe~ rt , where the subscript indicates the average over the initial state. The parameters 
used are 71 = 72 = 1 and u = 0. For the number and thermal paths, the initial photon number is 
(")o = 3. The parameters for the initial squeezed state are r = 1.2, a = 1. For the case of (a), the 
change in (n) upon photodetection is negative, while for (b) and (c), it is positive. This reflects the 
sub-Poissonian photon number distributions in (a) and the super-Poissonian distributions in (b) 
and (c). There also appears diffusive behavior in the no-count processes arising from homodyne 
detection in all of the three cases. The average behavior of damping (dashed curve) is consistent 
with sample paths. 

where p(t) := ^(1 — e~ Tt ) and we have used (78) in deriving the last equality. We note that 
Eq. (82) is a geometiric distribution. 

The Monte Carlo paths of (n) is shown in Fig. 2 (b). The average behavior showing 
an exponential damping in time is the same as the number state, while the change in the 
average photon number upon a photodetection is positive, reflecting the fact that the photon 
number distribution is super-Poissonian (see Eq. (21)) [16]. 

D. Squeezed state 

Finally, we consider the case in which a squeezed state is taken as the initial condition: 

\a,r) := D(a)S(r)\0), (83) 

where 

D(a) := e aaWa , (84) 
S(r) := e §( a2 - at2 ), (85) 
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with a being an arbitrary complex number and r an arbitrary real number which is called 
a squeezing parameter. 

To evaluate the generating function(al), we need to calculate 



(a, r\ : exp[Ka + K*a' + ua'a] : \a, r). 



(86) 



This is done in Appendix A with the result 

1 



(l + f(l 

x exp 



-2r 
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1 + e 2r 



(87) 



Thus, the generating functional and the generating function are given by 
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(88) 
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■ImKt 2Ima 
+ 



l + v t l + e 



2r 



2(Ima) 2 
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(89) 

l + e 2r 1 + I/t 

We may use these results to describe the joint photocount and homodyne measurement 
process. 

The Monte Carlo paths for initial squeezed state is shown in Fig. 2 (c). We take the 
initial parameters with r = 1.2 and a = 1. The change in the average photon number upon 
photodetection is negative because of the super-Poissonian photon-number distribution. 
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V. CONCLUSION 



We have introduced the measurement operators for the simultaneous measurement pro- 
cess of photoncounting and homodyne detection and derived the corresponding stochastic 
Schrodinger equation. These stochastic equations describe the time evolution of the quan- 
tum state under given measurement outcomes. The analytical expression of the conditional 
wave function is obtained and, using this expression, we have derived the probability den- 
sity function and generating functional of measurement records as explicit functions of the 
initial state of the system. We have also derived the marginal distribution functions of two 
output channels and generating functionals. We have applied these general results to four 
typical initial conditions: coherent, number, thermal, and squeezed states. For each of these 
initial states, we have obtained analytical expressions of the generating functional of the 
measurement records. For the initial coherent state, the two output records for the photon 
counting and homodyne channels are mutually independent, while for the other states the 
two outputs are statistically dependent. In this sense, this measurement process extracts 
non-classicality of the system in the form of the correlation of two output channels. We have 
performed Monte Carlo simulations of the average photon number, which show properties 
of both photon counting and homodyne detection, implying the particle-wave duality of the 
photon field. 
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Appendix A: Derivation of Eq. (87) 

In this appendix, we derive Eq. (87) based on the Q-function technique. 
The Q-function of state p is defined by [22] 




(Al) 



7T 
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This function is convenient for calculating the expectation values of anti-normally ordered 
operators. Let /(/?, /?*) be an arbitrary function of complex variables ft, (3* with an expression 

m?) = E (A2) 

r.s 

The anti-normally ordered operator of f((3,(3*) is defined by 

^[f(d,a^]:=J2frsa r ^y. (A3) 

Then, from the overcompleteness relation of coherent states, we have [22] 
tr(p^[/(a,at)]) = Yj rs Tr(pa r [ 



r. s 

2 



r,s 

d 2 pf(p,nQW,n- (A4) 

Equation (A4) implies that Q-function can be interpreted as a quasi-probability distribution 
function for anti-normally ordered operators. To exploit this property, we will calculate the 
Q-function of the squeezed state \oc,r) and the anti-normally ordered expression of the 
nomally-ordered operator : exp[nd + n*d) + ud^d] :. 

1. Q-function of the squeezed state 

To evaluate the Q-function of the squeezed state \at,r), we first prove following formula: 



S(r) = \l - e -^2^tanhr e -atalncoshr e i-tanhr_ 

coshr 



To show this, we differentiate the rhs with respect to r: 

— (rhs) 
dr 



t anh r {&^_ tanh r _ & t & ln cosh r s£ tank r ifl) 2 M£ tanh r „ - a+ a In cosh r „ ^- tanh i 



coshr \ 2 2 cosh 2 r 

_ e _(^ tanh ^_ a t atanhr ) e -atalncoshr e Vtanhr + tanh r^aU ln cosh r I \ ^tanhr 

\ 2 cosh r / 

a 2 — (a 1 ") 2 , , . ... 
^-^(rhs) (A6) 
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In deriving the last equality, we used the following relations: 



-7(a + ) 2 (n -^ t og e Aata e 7(at) 2 _ p \ 



e A (a + 27a 



(AT) 
(A8) 



Equation (A6) shows that the rhs satisfies the same differential equation for the lhs. Since 
(lhs) = (rhs) = I when r = 0, Eq. (A5) holds for arbitrary r. 
Using Eq. (A5) and 

D(-P)D{a) = e^*- a *^D(a - /3), 

we obtain 
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Therefore, the Q-function for the squeezed state is given by 



Q({3,n = -mD(a)S(r)\0)f 

7T 



7r cosh r 



exp 
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2. Ant i- normally ordered expression 



We derive the anti-normally-ordered expression of the operator 
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we have 



e utfa ; = (1 + j/fi 



Using 
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we obtain 
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3. Evaluation of Eq. (87) 



We can then evaluate Eq. (87). The relevant part is 



(a, r\£/[exp[K,'a + k'*ct + z/aa']]|a;, r) 
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(A16) 



Going back to original k and we obtain Eq. (87). 
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